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We present a minimal model for cuprate superconductors. At the unrestricted mean-field level, 
the model produces homogeneous superconductivity at large doping, striped superconductivity in 
the underdoped regime and various antiferromagnetic phases at low doping and for high tempera- 
tures. On the underdoped side, the superconductor is intrinsically inhomogeneous and global phase 
coherence is achieved through Josephson-like coupling of the superconducting stripes. The model is 
applied to calculate experimentally measurable ARPES spectra. 
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I. INTRODUCTION 

We are witnessing an increase of experimental evidence 
indicating that (charge and magnetic) incommensura- 
tion characterize the low-energy physics of underdoped 
cuprate superconductors, both above and below the crit- 
ical superconducting temperature Tc |l)]. This poses a 
challenging problem to theorists since these compounds 
appear to be at the verge of a multitude of different quan- 
tum ordered states that can be tuned by varying phys- 
ical parameters of the system. Theorists like to use the 
word "quantum criticality" to refer to this phenomenon. 
The truth is, however, that so far there is no rigorous 
theoretical framework that can explain unambiguously 
this variety of complex phenomena (superconductivity, 
magnetism, incommensuration, etc.), characterized by 
intrinsic nonlinearities producing large-scale sensitivities 
to small perturbations. 

In this paper we present a minimal model of high-Tc su- 
perconductors that clearly displays a variety of commen- 
surate and incommensurate competing thermodynamic 
phases. The advantage of our approach is that it is sim- 
ple, not subject to vague argumentation, and it allows 
one to rigorously and exhaustively explore a variety of 
physical observables. 

Crucial to the experimental findings have been neutron 
scattering techniques which probe the spin dynamics of 
the high- Tc compounds and suggest that different fam- 
ilies of cuprate superconductors share inhomogeneously 
spin and charge textured phases as their quantum states 
y. It is well-known that the stoichiometric (half- filled) 
compounds are antiferromagnetic (AF) Mott insulators 
as a result of strong electron interactions and it is upon 
charge doping that they display incommensuration. In- 
deed, in a recent paper [EJ we have presented a uni- 
fied theory for the commensurate resonance peak and 
low-energy incommensurate response observed in neu- 
tron scattering experiments. We ascribe both features 
to be purely magnetic in origin: They represent univer- 
sal features signaling the existence of an incommensu- 
rate spin state both below and above the superconduct- 



ing transition temperature. Our interpretation indicates 
that superconductivity is not the reason for the resonance 
peak, and that the incommensurate quantum state pro- 
vides a reference state for the underdoped cuprates. 

In previous work |g|,y we introduced two classes of mi- 
croscopically inhomogeneous models which captured the 
magnetic and pairing properties of underdoped cuprates. 
Starting from a generalized t-J model Hamiltonian in 
which appropriate terms mimic stripes, we found that 
inhomogeneous interactions that locally break magnetic 
SU(2) symmetry can induce substantial pair binding of 
holes in the thermodynamic limit. We showed that these 
models qualitatively reproduce the ARPES and neutron 
scattering data seen experimentally. Moreover, based on 
the phenomenology of our microscopic model we devel- 
oped a mean-field ( " Josephson spaghetti" ) model which 
provides a scenario for the macroscopic superconducting 
state. From our model Hamiltonian of random stripe sep- 
aration r and associated inter- and intra-stripe random 
Josephson coupling J(r) ^ \/r we obtained the exper- 
imentally observed relation Tc{x) ~ (Ji'r)) ex: [(r)]"^ — 
6{x), where x represents charge doping and S is the in- 
verse of a characteristic length scale associated with the 
incommensuration. 

Our previous numerical simulations have helped to elu- 
cidate a certain fraction of the underdoped cuprate puz- 
zle. Here we assume a different strategy complementary 
to the previous approach: We propose a minimal homo- 
geneous model based on the one-band repulsive Hubbard 
Hamiltonian on a square lattice. The attractive particle- 
particle singlet channel is included through the nearest 
neighbor attraction V [g| , which produces predominantly 
d-wave pairing close to half-filling. We solve this model 
at the mean-field level allowing all physical quantities to 
vary from one lattice site to another. In this way, quan- 
tum morphologies characterized by a certain correlation 
length ^ will appear when the minimum length of our 
supercell is larger than <^. 

The basic question we address in this paper is whether 
antiferromagnetic striped ordering and d-wave super- 
conductivity can coexist in a certain parameter range of 
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our model. In our opinion, this question is crucial for 
the understanding of the superconducting state in the 
underdoped materials. 

In the next Section we introduce the model and briefly 
describe the way we solve it. We then summarize the 
resulting competing quantum states in a phase diagram 
and discuss the spectral density. At the end of the paper, 
we review our main findings. 

Some of the results presented in this paper have been 
introduced in Ref. H ; the analysis of the spectral density 
is described here for the first time. 



II. MINIMAL MODEL 
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We consider here a minimal model that brings together 
stripes and superconductivity. The model is the two- 
dimensional one-band Hubbard Hamiltonian with an on- 
site repulsion U f^. Pairing correlations are introduced 
by including the nearest neighbor attraction V ||9[]. The 
effective minimal Hamiltonian is thus 



Ht-t'-u = -Y^ Ujc\„c^„ + ^Y ni^ri.i (1) 

H = Ht^t'^u + vYn,nj, (2) 



where the operator c\^ {'^ja) creates (annihilates) an 
electron with spin a on the lattice site i, and Ui = 
chc^-f + f^i\Ci\ represents the electron density on site i. 
The hopping iy equals t for nearest neighbors and t' for 
for the second-nearest neighbor sites i and j. For our 
computations, we use the unrestricted mean-field approx- 
imation to this Hamiltonian, 



Hmf = - XI *y"Ci<TS-a + U^rii^in^i) + {n^^)n^i 

+ IZ^UST^^J+H-C- , (3) 



where A^j = V{ciiCj-\) is the MF superconducting order 
parameter. The effect of V in our model is limited to the 
generation of superconducting correlations. We do not 
explicitly address the important issue of the microscopic 
origin of the attraction V ||] . 





FIG. 1. Typical example of density and superconducting 
order parameter profiles in a stripe state (here, period 17). 
The top two bar charts represent the site-dependent spin 
and charge densities, respectively. The contour plots indi- 
cate the sites with low (blue) and high (red) values of the 
corresponding densities. The bottom four plots show the 
values of the superconducting order parameters, defined as 
Aj " = (Ai,,.ight+Ai,!e/tTAi,up=FAi,doii.n)/4 for d-wave (ex- 
tended s-wave) order parameter on site i {U = 4t, V = — 0.9i, 
i' = 0). Different choices of parameters lead to qualitatively 
similar patterns, with stronger U leading to a stronger AF 
order and more attractive V causing the superconducting 
stripes to become wider and larger in amplitude. The doping 
level is 5.9%. 



A typical zero-temperature MF inhomogeneous solu- 
tion is shown in Fig. nl In the lowest energy configu- 
ration, the spin density develops a soliton-like AF anti- 
phase domain boundary — a stripe — at which the AF 
order parameter changes sign. At the domain boundary, 
the electronic charge density is depleted. The width of 
the domain wall, S,dW: decreases with increasing on-site 
repulsion U . However, for values of U that are not much 
larger than the hopping t, the charge per unit length of 
the optimal (the lowest energy) stripe remains the same 
and is close to unity near half- filling for t' = 0. The bond- 
centered stripes are favored relative to the site-centered 
ones, although the energy difference in our case is small 
due to the smooth charge distribution. 

Stripe formation is the result of the competition be- 
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tween antiferromagnetism (which can lead to charge con- 
finement) and delocahzation (driven by kinetic energy). 
Non-linear feedback is responsible for these complex pat- 
terns. A half-filled (one electron per site) antiferromag- 
net is the state with the lowest energy per electron. Upon 
small doping, the energy can be optimized by segregating 
the excess electrons or holes into domain walls or charged 
clusters, and keeping the bulk antiferromagnetism unper- 
turbed. The anti-phase domain walls are favored over the 
simple charged stripes as they allow the charge carriers 
to optimize their transverse fluctuations in the direction 
across stripes, thereby lowering their kinetic energy. The 
linear filling of the emergent stripes vary depending upon 
the specific nature of the band structure. 




FIG. 2. Linear filling of an isolated vertical stripe as a func- 
tion of next- nearest neighbor hopping, t' . Here [/ = 4t and 
V = 0. 



For different band structures the exact relation be- 
tween the doping x and inter-stripe distance, L{x), may 
change; however, any model whose ground state is AF 
at zero doping, can be expected to have AF stripes for 
a finite doping, with incommensuration proportional to 
the doping, l/L{x) ex x, near half-filling. For example, 
negative next-nearest neighbor hopping t' (relevant in the 
hole-doped cuprates [|lO[), modifies the stripe filling with- 
out compromising the stripe phase stability relative to 
commensurate AF at the MF level [|l| . The stripe filling 
is a monotonically decreasing function of the magnitude 
of t', with the filling 1/2 occurring when t' — — 0.35t 
(Figure 0). While filling-one stripes correspond to a cor- 
related insulator, the fractional filling stripes in the t-t'- 
U model (Eq. (Q)) are metallic, which can be understood 
in terms of the partial occupancy of the mid-gap band 
formed due to the stripes m^. In the case of insulating 
stripes (t' = 0) a threshold value of attraction V should 
be exceeded to generate superconductivity and hence to 
overcome the insulating gap ||l2]; however, in the case 
of metallic stripes one would expect that any attraction 



would yield superconductivity through the Cooper insta- 
bility ||l^. Indeed, this is what we find 0|. On the con- 
trary, the diagonal stripes which can also be the ground 
states of the Hubbard model, particularly at low dop- 
ings, always have a filling of one electron per Cu site, 
and hence are insulators. This makes diagonal stripes 
antagonistic to superconductivity pi[ , and also agrees 
with the experimental observations ]14| . 

From Figure |] it is clear that the superconducting or- 
der parameter A^ is maximized on the stripes and is 
not smooth (even within the stripe) due to the presence of 
the AF background. In addition to the dominant d-wave 
component, there is a small extended s-wave (s*) com- 
ponent generated on the stripe, which can be interpreted 
as a distortion of the d-wave at the level of about 10%. 
This happens because the symmetry of the lattice has 
been spontaneously broken by the stripes. For dopings 
less than about 10% (corresponding to L{x) > 10 lattice 
sites) the stripes have negligible overlap . In this regime, 
the amplitude of the superconducting order parameter on 
the stripes no longer depends upon the stripe-stripe sep- 
aration. For higher doping levels, an overlap between the 
superconducting order parameters on adjacent stripes is 
established, and for even higher doping the stripes "melt" 
and superconductivity becomes homogeneous, of a clas- 
sical BCS type 111. 



A. Phase diagram 




Doping, X 

FIG. 3. Schematic phase diagram obtained by super- 
imposing the antiferromagnetic (AF) / striped (ICAF) 
and the d-wave superconducting (SC) phase diagrams. 
In the intersection region we distinguish the subregions 
of Josephson-coupled striped superconductor (SSC), and 
non-superconducting "strange metal" (SM), which is neither 
a superconductor, nor a simple insulator. The upper bound- 
ary of the AF/ICAF corresponds to the weak pseudogap 
crossover, and the line between the pure AF/ICAF and the 
SM marks the strong pseudogap crossover. A detailed fi- 
nite-temperature study is required to precisely locate the left 
boundary of the SSC region, and hence to determine the order 
of the critical point P. 
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From our zero-temperature analysis of the coexistence 
of AF stripes (ICAF) and superconductivity , a sim- 
ple qualitative thermodynamic phase diagram emerges. 
In the conjectured phase diagram, we utilize the finite- 
temperature AF/ICAF phase diagram of the Hubbard 
model together with the superconducting (SC) phase di- 
agram of the t-V model. The SC phase diagram is ob- 
tained in the homogeneous MF ||], while the AF/ICAF 
phase boundary is constructed under the assumption of 
the second order phase transition between the homoge- 
neous and inhomogeneous states S . For a suitable choice 
of parameters, for instance U — 2t and V — — t, the SC 
and the AF/ICAF regions in the phase diagram inter- 
sect, as shown in Fig. y. The energy scale associated 
with the AF/ICAF region of the phase diagram is much 
larger than that of the SC part. Thus, one expects that 
only the SC phase boundary is modified when it passes 
through the AF/ICAF region. The central result of our 
work is that the superconductivity does not disappear in 
the region of the AF stripes, but rather becomes striped, 
with anisotropic superfluid stiffness. 

Based on familiar Josephson coupling physics, in the 
region of coexistence of superconductivity and stripes, 
we can expect a part that is a globally coherent striped 
superconductor (SSC). The rest of the intersection re- 
gion is covered by an exotic phase which, if it were per- 
fectly orientationally ordered, would be a superconduc- 
tor in one direction and a strongly-correlated insulator 
in the other. In reality, due to the meandering of the 
stripes and their break-up into finite segments |15[, the 
state is likely to be highly inhomogeneous and neither an 
insulator, nor a superconductor, but also not a simple 
metal. In agreement with the experimental attribution, 
we refer to this region as a "strange metal" (SM). The 
line separating the SM from the rest of the AF/ICAF 
region, in the context of the experiments, can be associ- 
ated with the crossover to the strong pseudogap regime, 
and corresponds to the opening of the local supercon- 
ducting gap. The high-temperature boundary separating 
AF/ICAF phases from the homogeneous state, marks the 
onset of the weak pseudogap. For small dopings, there 
is also a possibility of a transition from the vertical to 
diagonal stripes [0 . 



B. Spectral Density 

In this Section we demonstrate how our model can 
be applied to compute the energy spectrum of the sys- 
tem, which can then be compared with the experimental 
data |16|. At any temperature, the MF solution yields a 
self-consistent spectrum {£„} of Bogoliubov quasiparti- 
cles, which diagonalize the MF Hamiltonian in Eq. (O). 
Knowing how the electron operators are related to the 
Bogoliubov quasiparticles, one can compute the elec- 
tronic spectral densities for positive (particle) and nega- 



tive (hole) biases. Experimentally, angle resolved photoe- 
mission (ARPES) measures the electronic spectral den- 
sity integrated in the window of ±Ae around the Fermi 
energy. 

A e = 0.2 t 





AE = 0.3t 





Ae = 0.5t 





A e = 0.5 1 (non-interacting) 





FIG. 4. Computed ARPES spectra for various energy in- 
tegration windows, Ae. The Brillouin zone is defined as 
(0, 2tt) X (0, 27r), with the horizontal direction being along the 
stripes and the vertical direction perpendicular to the stripes. 
Notice that for small Ae the spectral weight is concentrated 
around (0,'7r), with the Fermi surface being gradually recon- 
structed with the increasing window of integration. The pa- 
rameters are U = 4t, V — —t, t' = — 0.2i, and doping 8.3%. 
The system size is 16 x 16, with two collinear stripes. For 
comparison we also show the case of free electrons with the 
same non-interacting band structure and doping. 



In Figure ^ we show the computed ARPES spectra for 
various energy integration windows. The symmetry of 
the spectrum is spontaneously broken due to the pres- 
ence of stripes (the stripes run along the horizontal k 
direction). For small integration window near the Fermi 
surface, the spectral weight is concentrated around the 
(0, tt) point. The reason is that the stripes gap out the 
flat parts ("diagonals") of the Fermi surface, while keep- 
ing the quasiparticles around (0, tt) gapless if t' < or 
weakly gapped if f' — 0. Since these are the quasipar- 
ticles that are primarily responsible for the formation of 
dr^2_y2 superconductivity, it is this particular structure of 
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the "stripe gap" that allows for the peaceful coexistence 
of stripes and superconductivity in our model. Indeed, 
in the absence of superconductivity {V = 0), the spectral 
patterns remain essentially unchanged, except for the en- 
hanced weight around (0,7r). 

For larger energy integration windows, the Fermi sur- 
face gradually "reconstructs," with the energy states 
around the diagonal reappearing when Ae exceeds the 
"stripe gap." 

Notice that due to twinning and the expected pres- 
ence of stripe domains in the real experimental systems 
the computed spectra have to be symmetrized. Similar 
results for the striped spectra (but without superconduc- 
tivity) have been obtained previously ||l^,|ll|,Q . Similarly, 
one can calculate temperature-dependent specific heat, 
entropy, spin susceptibility, among other experimentally 
measurable quantities, as we will report elsewhere. 



III. CONCLUSIONS 



In summary, we have presented a minimal model sup- 
porting the coexistence of incommensurate antiferromag- 
netism ( "stripe" order) and global anisotropic supercon- 
ductivity. Contrary to the common belief, these two 
order parameters can coexist and our calculation is a 
faithful realization of such a physical situation. At the 
same time, the stripe order provides a natural competing 
order parameter limiting the increase of the supercon- 
ducting transition temperature on the underdoped side 
of cuprates. 

Our model displays a variety of other competing ho- 
mogeneous and inhomogeneous thermodynamic phases. 
Based on the model, we constructed a phase diagram that 
captures many features of the superconducting cuprates. 
Finally, we computed the photoemission spectra which 
give a clear interpretation of experimental data. 

This work was supported by the U.S. DOE. 
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